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Abstract 

Natural disasters or attacks may disrupt infrastructure networks on a vast scale. Parts of the damaged 
network are interdependent, making it difficult to plan and optimally execute the recovery operations. To 
study how interdependencies affect the recovery schedule, we introduce a new discrete optimization problem 
where the goal is to minimize the total cost of installing (or recovering) a given network. This cost is 
I determined by the structure of the network and the sequence in which the nodes are installed. Namely, the 

^ cost of installing a node is a function of the number of its neighbors that have been installed before it. We 

' ' analyze the natural case where the cost function is decreasing and convex, and provide bounds on the cost 

of the optimal solution. We also show that all sequences have the same cost when the cost function is linear 
" ' and provide an upper bound on the cost of a random solution for an Erdos-Renyi random graph. Examining 

,— I the computational complexity, we show that the problem is NP-hard when the cost function is arbitrary. 

Finally, we provide a formulation as an integer program, an exact dynamic programming algorithm, and a 
greedy heuristic which gives high quality solutions. 
• Keywords: Infrastructure Networks; Disaster Recovery; Permutation Optimization; Linear Ordering Prob- 

''pi lem; Neighbor Aided Network Installation Problem. 

B 

1 Introduction 

Many vital infrastructure systems can be represented as networks, including transport, commu- 
nication and power networks. Large parts of those networks might become inoperable following 
C natural disasters, such as the Japanese 2011 earthquake near Sendai, to the extent that the cost of 

the recovery might be measured in hundreds of billions of dollars. Like earthquakes, wars or even 
terrorist attacks can cause massive network disruption. 

When faced with vast devastation, authorities must develop a plan or schedule for restoring 
the network. A particularly difficult and costly stage in the recovery is the initial stage because no 
J> infrastructure is available to support recovery operations. As the recovery progresses, previously 

installed nodes provide resources or bases that help reduce the cost of rebuilding their neighbors 
- a phenomenon that we call neighbor aid. For simplicity of analysis, we propose that during the 
reconstruction process all of the nodes and edges of the original network are to be re-installed. We 
assume that the cost of the edges is not determined by the schedule of the restoration. We examine 
the following question: How could the recovery schedule be optimized in order to reduce the total 
cost? 
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The phenomenon of neighbor aid apphes to many decision problems beyond disaster recovery. In 
fact it originally appeared in educational software design when one of the authors sought the optimal 
order in which to teach a list of vocabulary words. A word like neologism, for example, would be 
easier for students to learn if they already knew the etymologically adjacent words neophyte and 
epilogue. This means that the sequence with which words are taught could be optimized to make 
them easier to learn. 

Naturally, the problem we will formulate is related to scheduling problems such as single pro- 
cessor scheduling [2j , the linear ordering problem [Jj , the quadratic assignment problem [6] and the 
traveling salesman problem (TSP) [7J. Like TSP our problem asks for a schedule based on a graph 
structure, but unlike with TSP, the cost associated with visiting a given node could depend on all 
of the nodes visited before the given node. Our problem offers a new model for disaster recovery 
of networks. For other models see e.g. [U |3l [TOl |8]. 

2 General Formulation and Basic Results 

In this paper we introduce a description of the neighbor-aid phenomenon as a discrete optimization 
problem, which we call the Neighbor Aided Network Installation Problem (NANIP). An instance of 
NANIP is specified by a network G = (V, E) and a real- valued function / : No — )• The domain 
of / are non-negative integers up to the degree of the highest degree node in G. The goal is to find 
a permutation of the nodes that minimizes the total cost of the network installation. The cost of 
installing node vt £ V under a permutation o" of y is given by 

f{rivt,G,a)), 

where r{vt, G, a) is the number of nodes adjacent to vt in G that appear before vt in the permutation 
a. The total cost of installing G according to the permutation a is given by 

n 

GG{a) = Y,mvt,G,a)). (1) 
t=i 

The problem is illustrated in Fig. [T| 

NANIP could also be expressed as optimization over n-by-n permutation matrices 11 where A 
is the adjacency matrix of G: 

"li^E/ E^n^^n),, . (2) 

t=i \i=i / 

The inner sum computes the number of nodes adjacent to node t that are installed before t. 
Throughout the paper we use the standard notation n := \V\ and m := We assume that G is 
connected and undirected, unless we note otherwise. 

We begin with a preliminary lemma which establishes that all the arguments used in calculating 
the node costs must sum to m, the number of edges in the network. 

Lemma 1. For any network G, and any permutation a of the nodes of G, 

n 

J2f{vt,G,(j) = m. (3) 
t=i 
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Fig. 1: An illustration of the problem on a simple instance. When f{k) = j:^, nodes a and b are 
best installed early because they can reduce the cost of many of their neighbors. An optimal 
solution is the sequence a = (a, b, c, e, d, /, g, h) with cost 46 = 12 + 6 + 4 + 4 + 6 + 4 + 6 + 4. 

Proof. For every edge in G and any permutation cr, one of the two endpoints of the edge is the 
latter of the two to be installed under a. By definition, r{vt,G,a) is the number of edges whose 
latter endpoint is installed at time t. Thus summing up r{vt, G, a) for all t = 1, 2, . . . , n accounts 
for all the edges in G. □ 

One application of this lemma is the case of a linear cost function f{k) = ak + b, for some real 
numbers a and b. With such a function the optimization problem is trivial in that all installation 
permutations have the same cost. 

Corollary 2. For any G, when f[k) = ak+b for some real numbers a and b, then Gcio') = am+bn. 
Proof. 

n n n 

Cg(o") = ^ f{r{vt, G, cr)) = ^ {ar{vt, G,a) + b^ = " ''(^*' ^' cr) + bn = am + bn. 
t=i t=i t=i 

□ 

For some applications, it may be useful to know the cost of a sequence when the sequence 
is chosen at random. Our analysis of this case uses a frequently-used model of a random graph, 
namely the Erdos-Renyi model where any edge exists independently with probability p. 

Lemma 3. Let Gn,p be an Erdos-Renyi random graph on n nodes with the edge probability p. Let a 
be a permutation of the nodes ofGn,p chosen uniformly at random from the set of such permutations. 
For any real-valued function / : No — )• M+ and p > the expected cost satisfies 

^ n—l 
^ k=0 

Notice that the graph does not need to be connected in Lemma |3j 

Proof. Let a = (fi, . . . ,Vn) be a permutation on n nodes chosen uniformly at random from the 
set of all such permutations, where the label Vi is assigned to the i^^ node to appear in a. Then 
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when vt is installed at time t, the probability that exactly k adjacent nodes are already installed is 
— pY^^^^ ■ Hence the expected cost of installing vt is 



E 

A:=0 



t - 1 



p'{l-p)'-'-''f{k) 



for t = 1, . . . , n. Thus, the expected total cost of installing G according to a is 



/—I 7^— n V / i,_n +— 7. I 1 V / 



n— 1 n 



p'il-pY-'-'fik) 



t=l k=0 

n— 1 n— 1 



fc=0 t=k+l 



n— In— 1 ^ \ n— 1 n—1 / \ 

E E (I)f'(i - f r = E (T^/w E (I) (1 - f) 

A:=0 g=fc ^ ^ fc=0 ^ ''^^ q=k ^ ^ 



Since 



E 

t=0 



(1-x) 



for X G [0, 1) and any positive integer A;, we have that 



n-l ^ V oo ^ N 

E0u-.)'.E©a-.) 



g=fc ^ ^ <J=0 

Substituting this above and canceling terms we obtain the bound 

n-l 

p 



fc+1 



E 



fc=0 



□ 



An exact expression for '¥,\Cq^^(o)^ could be given using the 2-^1 hypergeometric function, 
where 

2Fi{a,b,c,z) = 2^ 



i=0 



Using Mathematica [S], we find that the expected total cost of installing Gn,p according to a is: 



n— 1 ^ n—1 ^ \ 



(l-p)^ 



k=0 



E 

fc=0 
n-l 



g=k 
(I-P)" 



P 



k+1 



n 



[l-pTi ^ ]2Fi{l,n + l,n + l-k,l-p) 



^ n—1 n—1 ^ \ 

- E /(^) - E /(^) ( - ^^r^'sFiCl' n + 1, n + 1 - A;, 1 - p) . 

^ fc=0 fc=0 ^ ^ 
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3 Decreasing Convex Cost Function 



The case of a decreasing convex cost function / is natural for many applications of NANIP because 
many times the neighbor aid phenomenon declines in significance as the number of neighbors 
increases. For a decreasing cost function / : No — )■ M+, we say / is decreasing convex if f{i) — 
f{i + 1) > /(j) — f{j + 1) for all j > i in Nq. (Notice that a differentiable, convex, non-increasing 
function f{x) : M M satisfies f{x) - f{x + \) > f{y) - f{y + 1) for all y > x.) 

It is illuminating to see the effect of convexity in / for a small graph, such as a 3-node path graph 
V = {^1,^2, vs} and E = {{vi, V2), (^2, ^^3)}- A sequence like {vi, V3, V2) has cost /(O) + /(O) + /(2) 
while a sequence like {vi, V2, v^) costs /(O) + /(I) + /(I). When / is decreasing convex, the second 
sequence is at least as cheap as the first one, since /(I) + /(I) < /(O) + /(2). 

For decreasing convex cost functions it is possible to find absolute lower bounds on the cost of 
installing a graph. To state this bound we extend the domain of / to non-integer values so that 
it is defined on all reals in [0,n — 1]. One approach that preserves convexity is to define f{q) for 
every positive non-integer q as the linear interpolation of f{[q\) and /([(?]). To derive the bound, 
we start with the case where only a subgraph H of the network needs to be installed. 

Theorem 4. Let H = {Vh,Eh) be a non-empty subgraph of G = {Vg,Eg)- Let a be any installa- 
tion sequence for G such that it first installs Vq \ Vh <ind then Vh ■ Suppose the nodes in Vq \ Vh 
have already been installed. Denote with Eqh the cut from Vq \ Vh to H : 

Egh = {{i,j} eEG:ieVG\VH,j eVn} . (4) 

Then for any decreasing convex cost function f, the total cost of installing Vh satisfies 

CH(a)>|yH|/f '^^' + '^^^' V (5) 



The scenario in Theorem |4] frequently arises in practical network recovery situations. The graph 
G may represent the total network of a region while the subgraph H represents just the damaged 
subgraph. As a result, the cost of rebuilding H is reduced by connections to existing infrastructure 
(the edges Egh)- 

Proof. The argument follows from Jensen's inequality 

f ( J2i=A \ ^ Yli=A fi^i) 

^yB-A + lJ - B-A + 1' ^ ' 

with Xi = r{vi, G,a), A = n + 1 — \ Vh\, and B = n. □ 

This bound suggests that the best installation sequence is one where the cost for every node is close 
to the average cost. A special case is where Vg \ Vh is a single node. Its installation cost must be 
/(O) and so we have a lower bound on the cost of G. 

Corollary 5. Let G = {V,E) be any graph and a any installation sequence. If f is non-increasing 
convex then 

CG(a) >/(0) + (n -!)/( — 
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An immediate application of this bound is to find the optimal solutions for tree graphs. On a 
tree T, the bound gives c(T) > /(O) + (n — 1)/(1). Consider any algorithm that picks a random 
starting node and then proceeds to install neighbors of installed nodes until all the nodes are 
installed. (An explicit example is stated later in the paper as Algorithm [2]) Such an algorithm 
would install T at a cost /(O) + (n — 1)/(1), which must be optimal by Corollary [s} 

We now proceed to give another bound by considering the degrees of the nodes. Assume that 
the nodes are labeled in order of increasing degree in G. Thus, if di gives the degree of node Vi, 
then di < d2 < ■ ■ ■ < dn- We consider the following relaxation of NANIP Q, in which we do not 
seek an installation sequence from which the values of r(vt,G,a) are deduced, but rather we seek 
a collection of numbers, pi, . . . ,pn, which satisfy constraints that are satisfied also by r{vt, G, a) in 
NANIP. More precisely for every i we require that pi is less than or equal to di, and we require that 
values of X^ILi^** ~ (compare to Lemma [T]). Thus we have the following minimization problem. 

n 

P = mmY,f{Pi) (7) 

i=l 
n 

S.t. ^^Pi = m 

i=l 

^ Pi 1^ di and E No for i = 1, 2, . . . , n. 



o) 6"' 
t 5-^ 

T3 3 
2 + 
1 



4 5 6 7 
node label 



Fig. 2: An illustration of the lower bound of ([9]). Gray areas indicate the degrees of the nodes {di) 
while the solid line shows an optimal solution (pj). Here s = 2 and r = 2. 

It is easy to find P exactly, giving an analytic lower bound on NANIP with decreasing convex 
/. Namely, P occurs when the first s of the nodes are installed with pi = di and the rest of the 
nodes are installed with pi equal to dg ox dg + 1, as follows. 

Theorem 6. For any G and decreasing convex f and any a we have the following hound. If 
di > m/n, let s = 0, otherwise let 



max Ik : {n — k)dk 

A;e{l,2,...,n} 



k 

E 



di < m 



(8) 



When s = 0, any optimal solution of ^ has n — (?7i — n[^J) elements in (pi)"=x equal to [^J, and 
m — n['^\ elements equal to [^J + 1; When s > 1, let r = m — (n — s)ds — X]i=i di, and an optimal 
solution is given by the following assignment: (i) pi = di for i = 1,2, ... ,s; (ii) for s < i < n: 
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Pi = ds + 1 for exactly r indices and pi = dg for exactly n — s — r indices (illustrated in Fig. 
Thus the cost of NANIP can he hounded from below hy 

C f^)>p=/(^-("^-"L^J))-^(LfJ) + ("^-^L^J)/(LfJ+l), whens = ^, 
\(n-r-s)/(4)+r/(4 + l) + n=iM)> whens>l. 

Proof. Convexity implies that for any integers a < 6 in Nq, 

/(a) + /(6)>/(a + l) + /(6-l). (10) 

where the equahty holds if and only if 6 — a = 1. In the case where di < m/n, ([s]) is well-defined 
and s > 1 because ndi < m and ^^11=1 — ^w-- Suppose that there is an optimal assignment 
{.Qi)^=i / (Pi)F=i- If there is j < s such that qj < dj then from ([s]) must exist / > s such that 



Qj < Qj' — 2. From (10) it follows that by assigning qj := qj + 1 and qj' := qj' — 1 we can improve 
the existing solution (or at worst, keep its cost unchanged). Hence, qj = dj for all j = 1, 2, . . . , s. 
Suppose now that there exists a pair of indexes s < k < k' < n such that \qk — Qk'\ ^ Applying 
the same transformation on q^ and qk' as in the previous case we can match or improve upon the 
existing solution. Hence for all s < £ < / < n we have \qi — qi'l = 1 • Finally, from the invariant 
Y17=i Qi = rn it follows that among the numbers Qs+i, qs+2, ■ ■ ■ ,Qn there are r of them equal to dg 
and n — s — r of them equal to ds + 1, which completes the proof. The case s = is simpler. From 
s = it follows that [^J < di < d2 < ■ ■ ■ < d^. Given a solution {pi)^^i^ if there were i,j such 
that \pi — Pj\ > 2 then we could match or improve {pi)f^i by transformations similar to above. □ 

As a corollary, one could show that the same bound applies in the case of non-strictly convex /, 
which might arise in applications. Given such an /, infinitesimally perturb it to obtain g which is 
strictly convex, where g{d) = f{d) + e*^, for some sufficiently small e > and d G {0, 1, . . . , n — 1}. 
The cost Ccify) with g is bounded by ([9]) and must be greater than the cost Cg{cf) with /, but the 
difference vanishes in the limit e — )• 0. 



4 Computation 

We now consider the hardness of solving NANIP and introduce three algorithms for the problem, 
an Integer Programming Formulation, a Dynamic Programming algorithm, and a fast heuristic. 

Theorem 7. The Neighbor Aided Network Installation Problem is NP-hard. 

Proof. NP-hardness can be established by reduction from the Max Independent Set problem. Given 
a graph G, we can find a maximal independent set by seeking the optimal solution to the instance 
of our problem with the same G and the function / defined by /(O) = and f{k) = 1 for all 
k ^ 0. Using that /, any optimal permutation found by solving our problem will necessarily have 
the maximum possible number of nodes with cost 0, and the rest with cost 1. The nodes with cost 
then necessarily form a (maximal) independent set, since if any two were adjacent then the latter 
of the two to be installed would have incurred a cost of /(I) = 1. □ 

Note that while this establishes that general NANIP is NP-hard, it does not apply to the case 
of a decreasing convex cost function, which remains an open problem. 
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4.1 Integer Programming Formulation 



NANIP with a decreasing convex cost / can be formulated as an integer program in the following 
way. For ah i £ V, let Xit = 1 iff a{t) = i, that is, Xa G {0, 1} denotes that node i is installed at time 
t G {1, . . . , n}. Observe that if and only if i is installed before j then for some T S {1, 2, . . . , n — 1}, 
Ylt=i ^it ~ Yld=i -^jt — 1- each edge {i, j) (considering direction) introduce a weight Eij G [0, 1] 
so that Eij + Eji = 1 for all i,j E V. We will ensure that at all optimal solutions, Eij takes integer 
values. Through Eij variables we express r(j, G, a), the number of neighbors of a node j installed 
before j: r{j,G,a) = J2ieN{j) ^ij- Lastly, let /^(x) define the downward slopping line through the 
points {d—1, f{d—l)) and {d, f{d)) with d G {1, 2, . . . , }. Collectively those lines define a polyhedron 
whose extreme points coincide with the values of f{d). Thus ld{x) = f{d) + {f{d) — f{d — 1)) (x—d). 
NANIP is the problem 



mm 

i=l 



Y,c^. (11) 



s.t. Cj>f{d) + {f{d)-f{d-l))\Y,Ei,-d\ j GV,d = l,2,...,n-l (12) 

.i&N{j) 



T 



Eij > i^jt -Xit) r = 1,2, . . . ,n - l,i G V,j G N{i) 

t=i 

Eij + Eji = 1 ieV,jeN{i) 

n n 

Y^^jt = 1> and Y^it = 1 ^ V,t = l,2,...,n 

t=l i=l 

Observe that because of the convexity of /, the points in the intersection of the inequalities 
(12) are exactly the points above the convex hull of {f{d)}^^^. At an optimal a, the value of Cj 
would minimize f{x) subject to x < YlieN{j) ^ij, so that the value of Cj = f{r{i, G, a)) . 

Using CPLEX version 12.2 (IBM Corp.) we were able to solve NANIP on tree graphs with 50 
nodes in less than one hour, and with longer running time as the number of edges increases. Con- 
vergence is greatly accelerated when the basic NANIP is generalized to allow for realistic variation 
in the cost of a node (/(d) becomes fi{d) for each i G V^). The speedup is likely due to the removal 
of symmetry in the problem, which is very high. For instance, in any graph on 2 or more nodes 
and any a, swapping the first and second nodes in a does not change its cost. 



4.2 Dynamic Programming Algorithm 

Consider a generalization of NANIP which we call the Subgraph- Aided Network Installation Prob- 
lem (SANIP). SANIP allows the cost r{vt, G, a) of a node at time t to be an arbitrary function of 
the subgraph induced by the first t — 1 nodes of a, rather than just the neighbors of the node. This 
cost may even depend upon the installation time because it is the number of nodes |y(-ff)|. 

SANIP instances can be solved by the following Dynamic Programming Algorithm, Algorithm[T| 
which exploits the fact that the cost of any node u depends only upon H and not how H was 
installed. Instead of considering the space of permutations on n nodes, which has size n!, it 
needs to consider merely the set of subsets of n nodes, which grows as 2". Concretely, to find an 
optimal installation of a subgraph on nodes T C V, it considers every decomposition of the form 
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T = S D {u}. The algorithm maintains a table TOPT indexed by subsets S <Z V , and stores a 
least-cost sequence for installing S along with its cost. The maximum size of the TOPT table is 
(^"2)- the pseudocode, the symbol V on sequences denotes concatenation. 

Algorithm 1 Dynamic Programming Algorithm for Construction of G 
TOPT[{u]] ^ (u) for all ueV. 
for all t e {2 . . . n} do 
TOPT' ^ 
for all S G TOPT do 
for all u G y \ S" do 
T ^ S\J{u] 

if T ^ TOPT' or cost{TOPT[S]) + c(u, S) < cost{TOPT'[T]) then 
TOPT'lT] <- TOPT[S] V (n) 
free TOPT 
TOPT ^ TOPT' 
Output TOPT[V] 



The algorithm uses 0(2**/^) space (the table size peaks at iteration n/2 and then declines) and 
0(n2") time. We found that on a desktop computer the computation is tractable up to n ~ 25. 
While the space requirement might be improved, the running time cannot be improved beyond 
i7(n2"). The lower bound occurs because it is possible to have an instance where the cost of 
installing some node u is arbitrarily small if u is sequenced immediately after some subgraph H 
of G. Observe also that the specification of SANIP would in general use 0(n2") space. To solve 
instances of larger size this algorithm could be replaced by a dynamical programming algorithm 
with a finite horizon. 

4.3 A Linear- Time Heuristic 

Exponential time algorithms can be avoided if suboptimal solutions are acceptable, which is often 
the case in practice. We will show that a simple greedy algorithm runs in linear time and we will 
show results of simulations comparing the results of the greedy algorithm with those of Algorithm [T} 
Algorithm [2] below preferentially installs nodes that have the lowest installation cost, which is to 
say, nodes that have the most installed neighbors. The cost of the solution found by this heuristic 
can be bounded when / is decreasing convex, as follows. In the worst case, the first node has cost 
/(O), and all other nodes are constructed at cost /(I). Often the bound can be made tighter, as 
follows. 

Algorithm 2 Greedy Algorithm for NANIP 

If n=0 Output (Taux 

Caux Random u £ V 
while I (Taux I < n do 

u G argmax^gyy^^^ {f{.r{u, G, daux)) : u £ V \ <Taux}, resolving ties arbitrarily. 

0"aux ^ faux V U. 
Output CJaux 
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Proposition 8. Let f be decreasing convex and G a graph with node degrees di 
and n > 3. The cost of the solution a found by Algorithmic is bounded as 



< do < ■ ■ ■ < dr, 



Ccicj) < /(O) + + fiq) + /(4+3) + • • • + f{dn) 



(13) 



where s G {0, 1, . . . , n — 1} is the largest possible, subject to the constraints (i) q £ {1,2, ... , ds+2} 
and (a) s + q + ds+3 + ■ ■ ■ + dn = m. 



In (13), the terms after s/(l) occur only if dn > 1. 

The proof, using techniques similar to the proof of Theorem [6| is elementary. Given that we 
can bound from below the cost of the optimal solution (Corollary [5]), one can even compute an 
optimality gap for any instance of NANIP. 

The practical performance of Algorithm [2] is illustrated by simulations in Fig.[3j Algorithm [2] is 
within 5% of the optimum in our simulations in networks of low, medium and high edge densities. 
Therefore, and because it runs in linear time. Algorithm [2] might be appropriate to large-scale 
practical instances of NANIP. Also shown in Fig. [3] are the degree upper bound and the relaxation 
bound. The simulations take a random connected graph on 15 nodes and m edges and the convex 
function f{i) = for i G Nq (the results are qualitatively similar for dozens of convex cost 
functions we examined.) For each value of m, 5 random connected graphs were generated by 
taking a random tree on n nodes and adding m — (n — 1) edges to it at random. The greedy 
algorithm was scored based on the average over 10 runs on each graph. 
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Fig. 3: The performance of bounds and of Algorithm [2] (Greedy) as a function of the number of 
edges, m. The greedy heuristic usually finds solutions indistinguishable from the optimum. 
The degree upper bound is from (13) and the relaxation bound is from 



In addition to this greedy heuristic, we have examined a heuristic which installs nodes based 
on their degree starting with the highest degree nodes. In our simulations, this heuristic did not 
outperform the cost-greedy heuristic (Alg. [2]) even when the cost function / was decreasing concave 
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(detailed results are not shown). Thus, Alg. [2] seems to be appropriate for the case of a decreasing 
concave /. 

5 Conclusion 

This paper introduces a new discrete optimization problem, the Neighbor- Aided Network Installa- 
tion Problem (NANIP). We found that NANIP and its generalization could be solved using integer 
and dynamic programming, while a fast greedy approach exists for the case of convex decreasing 
/. We therefore find that recovery operations of simple infrastructure networks could be planned 
by a simple rule: every step of the recovery should focus on the most accessible of the damaged 
network nodes. 
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